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Various Interpretations of the Stored and the 
Radiated Energy Density 

Miloslav Capek, Member, IEEE, and Lukas Jelinek 

Abstract —Three contradictory but state-of-the-art concepts for 
defining and evaluating stored electromagnetic energy are treated 
in this communication, and are collated with the widely accepted 
definition of stored energy, which is the total energy minus the 
radiated energy. All three concepts are compared, and the results 
are discussed on an example of a dominant spherical mode, which 
is known to yield dissimilar results for the concepts dealt with 
here. It is shown that various definitions of stored energy density 
immanently imply diverse meanings of the term “radiation”. 

Index Terms —Antenna theory, electromagnetic theory, electri¬ 
cally small antennas, Q factor. 

I. Introduction 

The evaluation of stored electromagnetic energy and its den¬ 
sity is one of the old but as yet unsolved problems of classical 
electromagnetism. This is true despite its straightforward and 
generally accepted definition; stored electromagnetic energy 
is that part of the total electromagnetic energy that is, in 
comparison with the radiated energy, bound to the sources 
of the field, being unable to escape towards infinity. In other 
words, the stored electromagnetic energy is the difference 
between the total and the radiated electromagnetic energies, 
representing a “rest mass” seen by the force exerted by the 
source. 

In the case of a static field and a quasi-static field, the 
evaluation of the stored energy is immediate, as it is just equal 
to the total energy m, the physical interpretation of which is 
directly inferred from Poynting’s theorem HI. The problem 
arises for fields generated by general radiators. One of the 
core problems is that within the time harmonic steady state, 
the total energy is infinite HI. This infinite energy is contained 
in the radiation field, and the goal of the evaluation of the 
stored energy is to subtract this infinite radiation energy. The 
second problem is that Poynting’s theorem gives us no clue 
of what radiation energy really is. One can only rely on some 
general properties like the positive semi-definiteness of the 
energy, and the fact that the far field of the radiator placed in 
the lossless media carries solely radiation energy El. 

Probably the first treatment of stored electromagnetic energy 
dates back to the work of Bateman El, who pointed out 
that the electromagnetic energy in a vacuum does not in 
general move with the speed of light in a vacuum, being 
slowed down by a kind of “rest mass”. This “rest mass” 
vanishes only in the case of pure radiation fields, i.e. only 
at an infinite distance from a finite source. This work has 
however been forgotten and it was not until recent years that 
it was recovered and generalized by Kaiser El, into the form 
of stored electromagnetic energy density. 
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In parallel, the problem of stored energy has also been 
extensively studied in the community of electrical engineering, 
mostly in connection with antennas. In particular, antenna 
designers commonly aim at the lowest energy storage in order 
to maximize the radiation efficiency of an antenna, and for 
this purpose one actually encounters the problem of stored 
energy evaluation El - To that point, Chu ID proposed a circuit 
equivalent of the spherical modes, and with its help subtracted 
the radiation energy and established fundamental lower bounds 
of the radiation quality factor. His method has been generalized 
by several works of Thai m, a, a. Radiation energy 
subtraction has also been attempted directly on the field level. 
Among the most prominent works we mention Ho), HU, HI, 
HI, HI, and also na, in which subtraction is taken from 
another point of view, claiming that separation of total energy 
into stored energy and radiated energy is unphysical, as it 
cannot be derived from Maxwell’s equations. 

The works of Bateman El, Kaiser El, and Collin HI 
attempted to define stored energy locally via its density, while 
the works of Rhodes ifTfil . Yaghjian and Best ifTTll . Vanden- 
bosch ns, and Gustafsson and Jonsson EH operated solely 
with stored energy as a whole. The local approach however 
offers certain benefits: first, radiation energy extraction is 
carried out at every space-time point, avoiding cumbersome 
operations with ill-defined infinite space integrations ns, and 
second, the definition via the density is also more physical, as 
within classical relativistic theory all laws should be formu¬ 
lated strictly locally. 

The concepts mentioned above yield similar values of stored 
energy for common radiators. However, there exist specific 
cases for which the methods are profoundly different. One 
such case is the dominant TE spherical mode, which has been 
shown to lead to apparently incorrect negative stored energies 
within some evaluation schemes na. 

This paper has two main purposes: to compare different 
interpretations of radiated energy, and to recall the possibility 
of a local definition of stored energy, which, despite its 
appealing properties, has not found its place within the antenna 
community. 

The paper is organized as follows. The necessary definitions 
and nomenclature are introduced in Sec. |II] Sec. |m] recalls 
several radiation energy extraction techniques. The techniques 
are compared on an example of the dominant spherical mode 
in Sec. HYi The results are discussed in Sec. m Conclusions 
are drawn in Sec. El 

H. Definitions 

In this paper, we will strictly omit dispersive media, 
in which the concept of stored energy is problematic, 
even without radiation ll20l . Furthermore, as is common 
in the theory of electromagnetic radiators, this paper will 
deal with time harmonic fields represented by field pha- 
sors F {r,uj) at angular frequency w, which relates to the 
time domain quantities as F (r,t) = Re{i^ (f, w) exp(jwf)}. 
We will also use cycle mean averages, which are defined 
as (/(t)) = (1/T) / (i) df, with T = 2TTlio, and are 

widely employed in the case of power quantities, where 
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{T (r, t) • Q (r, t)) = (1/2) Re {F (r,a;) • (r, a;)} within 

the time harmonic domain ED, and in which symbol * denotes 
complex conjugate. Note, however, that the results detailed 
below and summarized in Table |I] hold for a general time 
domain held. 

For the purpose of comparing various concepts of stored 
energy, we will advantageously use a dimensionless quality 
factor, dehned as 


_ OjWsto 
~ Pr 


( 1 ) 


where 

M4to = J {Wsto {r, t)) dV (2) 


is the total cycle mean stored energy, dehned via its density 
Wsto i'f'jt), where 


Pr = (n){£ {r, t) xTi (r, t)) ■ dS 


Re{£;(r,a;) x H* (r,w)}- dS 


(3) 


is the cycle mean radiated power m. 


III. Various Definitions of Stored Energy Density 
AND Radiated Energy Density 

This section briehy reviews three major dehnitions of stored 
electromagnetic energy density used in the literature. The 
three concepts are summarized in Table and their general 
properties are discussed in the following section. 


B. The Concept of Kaiser-Bateman 

A very interesting way of evaluating stored energy was 
proposed by Kaiser i), generalizing the previous work of 
Bateman 0. Within this concept, the relativistic energy- 
momentum relation m is used to dehne the stored energy 
density exposed in the second row of Table |I] 

To the best of the authors’ knowledge, this is the hrst time 
that stored energy density has been dehned strictly locally 
with no reference to the position of the sources, i.e. radiation 
energy extraction is carried out locally at every point of the 
space-time. 

Using time harmonic helds, the cycle mean stored energy 
within this concept can be written explicitly as 

(w^o (6) 

in which 


= ^(e\\E {r,uj)f + ^i\\H {r,uj) 
+ ^R.e|(eE; {r,u}) ■ E{r,w) 

+ pH {r,u}) ■ H 


— (Re {E {r, oj) x H* (r, w)} 


Co 


-f Re{E;(r,a;) x (r, w) ) 


(7) 


with the corresponding quality factor 


A. The Concept of Collin-Rothschild 

The classical scheme for radiation energy extraction was 
dehned by Collin and Rothschild in ifTOll . and was later rehned 
in HID into the form of energy density exposed in the hrst 
column of Table The corresponding cycle mean energy 
density for a time harmonic held can also easily be written 
as 

(rufR ^r,uj)) {e\\E (r,a.) f -f p\\H {r, cu) f) 

^ 1 (4) 

— -—Re {E (r, oj) x H* (r, w)} • no, 


where tiq is the normal to the far held wave-front. 

This radiation energy extraction is the most common method 
used in the literature 0, despite its immediate dehciency 
of using riQ as the direction of the power How. This poses 
no problem for specihc geometries (pure modes in separable 
coordinate systems), but cannot suffice in general. Any general 
radiator will clearly in its near held emit radiation in directions 
different from the far held. 

Eor the sake of comparison, the quantity 




Pr 


(5) 


is dehned as advantageously normalized total stored energy 
within this concept. 


KB 


qKB _ ^^sto 
^ Pr 


( 8 ) 


Two crucial differences immediately appear when compar¬ 
ing 0 and (|^. Eirst, in (|^ the entire power how is subtracted 
from the held energy, while in Q only the power how along 
the direction of wave-front at inhnity is subtracted. Second, 
in the subtraction in done is squares, while in 0 the 
subtraction is direct. As a result of the squared subtraction 
in the cycle mean cannot be simply performed a priori. 

With respect to this concept of Kaiser-Bateman, it is worth 
mentioning a situation in which the general time domain 
dehnition, see Table |Ij can be greatly simplihed 11. This 
happens in the case when S (r, t) ■ H (r, t) = 0, which is the 
case for the example in this paper, and it is not rare even in 
realistic situations (at least in an approximate sense). Under 
this specihc condition, the dehnition from the second row of 
Table U can be rewritten as 


w 


KB 

sto 


(r,t) 




(9) 


which is a rather curious form. Not being the absolute value, 
(|^ would correspond to the energy excess appearing in the 
complex Poynting’s theorem 0. The absolute value, however, 
makes it (according to Kaiser and Bateman) the stored energy 
density. 
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TABLE I 

Various concepts of extraction and subtraction of radiated energy. 


Concept 

Wato (t, i) 

Wrad (r, t) 

Collin and Rothschild 1101 

U {r,t) - Wrad {r,t) 

S{r,t) 

• no 

CO 

Kaiser 0, Bateman 0 

(r, t) - (r, t) 

S (r,t) 

CO 

Rhodes 1161. Yaghjian and Best 1171. Vandenbosch 1141. Gustafsson and Jonsson 

U {r,t) - Wrad {r,t) 


W (r, r) = |eo||f (»■, t) IP-1-(r, t) ||2, S {r,t) = £ (r,t) X'H (r,t), tF {r,t) = lim {r£{r,t+-f)) 

r—^oo ' ' 0 / / 


C. The Concept of Rhodes 

Another well-established scheme for stored energy evalua¬ 
tion was used by Rhodes M, was generalized by Yaghjian 
El and was later reworked into the source concept by Van- 
denbosh IIT4ll for radiators of arbitrary shape. The definition 
of stored energy density within this concept is exposed in the 
third row of Table The formula can also easily be rewritten 
for time harmonic fields as 

= i ||£; (r,a;)f + M ||Tf (r, 

( 10 ) 

in which F = lim (rE exp(jfcr)). The corresponding nor- 

r—¥oo 

malization is dehned as 




Fig. 1. Sketch of the coordinate system that is used throughout the paper 
(a) and a sketch of the TEio current on a spherical shell of radius a. (b) The 
input current is normalized to /q = 1A with respect to the light blue contour. 


( 11 ) 

-fr 


IV. Stored Energy and Its Density For the 
Dominant Spherical TE Mode 

The dominant spherical TE mode is defined as the field 
generated by the cuiTent density HI 


sin (i?) 

J = 2a ^ 


H = 


kah^j^\ka) f . ^ ^ 

gfc (^-Ji (ka) cos (t?)ro 


kr 

+ tTio(fcr)-ii(fcr) ^i,^(^),, 


and for r > a 


E = - ^kaii (ka) (kr) sin (i?) 

H ( — hp^ {ka) cos (t?) Tq 

kr V 

,(2) 


,(2) 


where k = uj/cq is the free space wave-number and cq is the 
speed of light. 


speed of light. 

The total energy of the fields 1E*°* (in the form of the__ 

sponding quality factor) can be evaluated in a straightforward 

m^innpr iic 


( 12 ) 


manner as 

with 


. ,TJ/tot 

/^tot _ I /^tot 

^ ~ p ~ ^int “k Vext 


(15) 


flowing on a spherical shell of radius a, see Fig. [T] and it pro¬ 
vides interesting testing grounds for stored energy evaluation 
ifTSl . The fields generated by this source read fT2\ for r < o 

E = — ^^ka h^^^ (fca) ji (kr) sin (i?) ipQ, (13a) 


t _ 

'^int — Tj 


ka 

= j Qint‘ (kr) dkr 


ka 


(ka) T , 

—TJf^ J [ikr)^\ji{kr)\^+^\ji{kr)\^ 

0 


(13b) 

-(- fcr jo {kr) - 



tot _ _ 7 

C/ext p J 

(14a) 

ka 

OO 


(16a) 


Qex‘t (kr) dkr 


kail{ka} 
kr 

+ (14b) 


1 

" 2 
-f 


OO 

J (^{kr)^ {kr) + 2 {kr) 

';a 

kr {kr) — {kr) ^ dkr, 


(16b) 












































JOURNAL OF LTbX CLASS FILES, VOL. XX, NO. XX, XXX 20XX 


4 



ka 10 



Fig. 2. The normalized total energy density and the normalized stored energy 
density of the spherical TEio mode for ka = 6. Several approaches to obtain 
the stored energy density are depicted. The energy jump is given by the 
presence of the current shell at r = a. The surface of the shaded area is 
equal to ka, which is in exact correspondence with tH- 


where subscript “int” denotes energies for r < a and subscript 
“ext” denotes the energies for r > a. The integration in ( |16a| i- 
(16b I can be carried out analytically (though it will lead to 


Qext cxd)’ but we rather leave out the possibility to study 
the radial energy density represented by Q (kr). 

Similarly, we can evaluate the stored energy and its ra¬ 
dial density within the scheme of Collin-Rothschild, Kaiser- 
Bateman and Rhodes as 

gCR/KB/Rh _ gCR/KB/Rh gCR/KB/Rh 


with 

ka 

gCR/KB/Rh ^ r j,CR/KB,Rh 

0 

oo 

^CR/KB/Rh ^ I gCR/KB/Rh 
ka 

The densities 2-°^ (kr), 2*°^ (kr), 

(fcr) are depicted for fca = 6 in Fig. while 
the and QCR/KB/Rh depicted in Fig. I as functions 
of ka. 


V. Discussion 

There are several important observations in Table 0 Fig- in 
and Fig. that will be discussed separately in the following 
subsections. 

A. Over-subtraction 

The curves in Fig. [^reveal that radiation energy subtraction 
inside the circumscribing sphere brings serious issues, the 
significance of which will grow with the electrical size of the 
radiator. 

Particularly, there is an observable difference between the 
stored energy density of Collin-Rothschild and that of Kaiser- 
Bateman. While the scheme of Collin-Rothschild does not 
subtract any radiation energy in the internal region (the power 


Fig. 3. The normalized total energy and the normalized stored energy as a 
function of ka. Several approaches to obtain the stored energy are depicted. 


flow in the radial direction is strictly zero), the scheme of 
Kaiser-Bateman subtracts the entire Poynting’s vector there, 
which is indeed non-zero also inside the sphere. A comparison 
of the stored energy density in Kaiser-Bateman’s scheme with 
the total energy density reveals that there is an energy loss 
in the internal region provided by angular components of the 
Poynting’s vector. This energy is however not radiated out of 
the sphere. This subtraction is clearly incorrect, as it will also 
be performed in the case of a spherical cavity, which does 
not contain radiation, at least not in the classical sense of an 
energy reaching infinity. 

A different form of over-subtraction also burdens the 
scheme of Rhodes, but some care should be taken with the dif¬ 
ferent interpretation of this scheme by various authors, namely, 
Rhodes ifT^ . Yaghjian and Best ifTTl , Vandenbosch ifT^ . and 
Gustafsson ifTSl . In particular, Rhodes ifT^ divides the entire 
space similarly as in this paper, see ([nil, thus obtaining the 
internal stored energy and the external stored energy. Radiation 
subtraction according to the third row of Table[I|is used only in 
the external region, making this scheme identical to the scheme 
of Collin-Rothschild ifTOl . By contrast, Gustafsson’s approach 
ifTSll uses subtraction everywhere, but adds ka to the resulting 
stored energy, which for spherical radiators gives the same 
result as the method of Collin-Rothschild, thought it generally 
differs. Finally, the approach of Yaghjian and Best ifTTll and 
the approach of Vandebosch IH use subtraction everywhere, 
with no further compensation. This leads to over-subtraction 
of the order ka in normalized scale that is used. 


B. Positive semi-definiteness 

Figure [^ and Fig. [^ show that the stored energy density and 
the total stored energy are both positively semi-definite within 
the scheme of Collin-Rothschild and Kaiser-Bateman. In fact, 
the scheme of Kaiser-Bateman is manifestly positively semi- 
definite H, and general positive semi-definiteness can also be 
expected from the scheme of Collin-Rothschild im. 

By contrast, the scheme of Rhodes within the paradigm of 
Yaghjian, Best and Vandebosch can lead to negative stored 
energies, which is clearly unphysical, see also lfT9l . This 
issue however comes as no surprise when we realize that 
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for a current distribution independent of frequency (the case 
presented in this paper), the aforementioned stored energy can 
be expressed as ll^ 

where X is the reactance seen by the sources of the field, 
where /q corresponds to the appropriate current normalization, 
and where dXjdijj is well known to reach negative values in 
radiating systems ll24l . The source of this negativity problem 
is clearly over-subtraction of the radiation energy, see Fig. 

C. Coordinate dependence 

One of the basic requirements on a valid physical quantity is 
its independence from an absolute coordinate system. Unfortu¬ 
nately, this is not satisfied in the case of the first and the third 
line of Table [I] where the radiation energy definition explicitly 
refers to absolute coordinates. A change of the coordinate 
origin then leads to a change in the value of the radiated and 
stored energy, which contradicts the idea of energy storage 
being a property of the radiator. This problem of the scheme 
of Rhodes is in fact well known, and some ways of minimizing 
it have already been published iflTl . ifTSll 

VI. Conclusion 

This communication has reviewed and discussed three up- 
to-date concepts of stored electromagnetic energy density. 
It has been shown on a particular example that, although 
sound, all three concepts yield results that to a certain degree 
contradict physical reality. It has been shown that the problems 
result from an improper definition of radiated energy density. 
Particularly, the concepts that have been discussed failed at 
least in some of the following prerequisites for a physically 
meaningful definition, which should 

• be strictly local, 

• be coordinate independent, 

• be gauge invariant, 

• give a positively semi-definite energy density, 

• give a zero value everywhere in a closed cavity. 

As a result, the correct definition of radiated energy density 
remains an open question, despite its long history and despite 
being a fundamental question in the classical theory of elec¬ 
trodynamics. 


[8] -, “Polarization, gain, and Q for small antennas,” IEEE Trans. 

Antennas Propag., vol. 59, no. 12, pp. 4844—4848, Dec. 2011. 

[9] -, “Q bounds for arbitrary small antennas: A circuit approach,” IEEE 

Trans. Antennas Propag., vol. 60, no. 7, pp. 3120-3128, July 2012. 

[10] R. E. Collin and S. Rothschild, “Evaluation of antenna Q,” IEEE Trans. 
Antennas Propag., vol. 12, no. 1, pp. 23-27, Jan. 1964. 

[11] R. E. Collin, “Minimum Q of small antennas,” Journal of Electromag¬ 
netic Waves and Applications, vol. 12, no. 10, pp. 1369-1393, 1998. 

[12] D. R. Rhodes, “On the stored energy of planar apertures,” IEEE Trans. 
Antennas Propag., vol. 14, no. 6, pp. 676-684, Nov. 1966. 

[13] -, “A reactance theorem,” Proc. R. Soc. Lond. A., vol. 353, pp. 1-10, 

Feb. 1977. 

[14] G. A. E. Vandenbosch, “Reactive energies, impedance, and Q factor of 
radiating structures,” IEEE Trans. Antennas Propag., vol. 58, no. 4, pp. 
1112-1127, Apr. 2010. 

[15] S. M. Mikki and Y. Antar, “A theory of antenna electromagnetic near 
field - pai't I,” IEEE Trans. Antennas Propag., vol. 59, no. 12, pp. 4691- 
4705, Dec. 2011. 

[16] D. R. Rhodes, “Observable stored energies of electromagnetic systems,” 
J. Franklin Inst., vol. 302, no. 3, pp. 225-237, 1976. 

[17] A. D. Yaghjian and S. R. Best, “Impedance, bandwidth and Q of 
antennas,” IEEE Trans. Antennas Propag., vol. 53, no. 4, pp. 1298- 
1324, April 2005. 

[18] M. Gustafsson and B. L. G. Jonsson, “Stored electromagnetic energy 
and antenna Q,” Prog. Electromagn. Res., vol. 150, pp. 13-27, 2014. 

[19] M. Gustafsson, M. Cismasu, and B. L. G. Jonsson, “Physical bounds and 
optimal cun'ents on antennas,” IEEE Trans. Antennas Propag., vol. 60, 
no. 6, pp. 2672-2681, June 2012. 

[20] L. D. Landau, E. M. Lifshitz, and L. P. Pitaevskii, Electrodynamics of 
Continuous Media, 2nd ed. Butterworth-Heinemann, 1979. 

[21] R. F. Harrington, Time-Harmonic Electromagnetic Fields, 2nd ed. John 
Wiley - IEEE Press, 2001. 

[22] J. A. Stratton, Electromagnetic Theory. John Wiley - IEEE Press, 2007. 

[23] M. Capek, L. Jelinek, P. Hazdra, and J. Eichler, “The measurable Q 
factor and observable energies of radiating structures,” IEEE Trans. 
Antennas Propag., vol. 62, no. 1, pp. 311-318, Jan. 2014. 

[24] S. R. Best, “The Foster reactance theorem and quality factor for 
antennas,” IEEE Antennas Wireless Propag. Lett., vol. 3, no. 1, pp. 306- 
309, Dec. 2004. 


References 

[1] J. D. Jackson, Classical Electrodynamics, 3rd ed. John Wiley, 1998. 

[2] J. Schwinger, L. L. DeRaad, K. A. Milton, and T. W.-y., Classical 
Electrodynamics. Westview Press, 1998. 

[3] H. Bateman, The mathematical analysis of electrical and optical wave- 
motion on the basis of Maxwell’s equations. Forgotten Books, 2010. 

[4] G. Kaiser, “Electromagnetic inertia, reactive energy and energy flow 
velocity,” J. Phys. A.: Math. Theor, vol. 44, pp. 1-15, 2011. 

[5] J. L. Volakis, C. Chen, and K. Fujimoto, Small Antennas: Miniaturiza¬ 
tion Techniques & Applications. McGraw-Hill, 2010. 

[6] L. J. Chu, “Physical limitations of omni-directional antennas,” J. Appl. 
Phys., vol. 19, pp. 1163-1175, 1948. 

[7] H. L. Thai, “New radiation Q limits for spherical wire antennas,” IEEE 
Trans. Antennas Propag., vol. 54, no. 10, pp. 2757-2763, Oct. 2006. 



